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Intracellular cargos which are transported by molecular motors move stochastically along cy- 
toskeleton filaments. In particular for bidirectionally transported cargos it is an open question 
whether the characteristics of their motion can result from pure stochastic fluctuations or whether 
some coordination of the motors is needed. The results of a mean-field model of cargo-motors dy¬ 
namics, which was proposed by Muller et a/. [I] suggest the existence of high motility states which 
would result from a stochastic tug-of-war. Here we analyze a non-mean field extension of their 
model, that takes explicitly the position of each motor into account. We find that high motility 
states then disappear. We consider also a mutual motor-motor activation, as an explicit mechanism 
of motor coordination. We show that the results of the mean-field model are recovered only in case 
of a strong motor-motor activation in the limit of a high number of motors. 


INTRODUCTION 

Almost every cellular function is related to transport 
processes. These are necessary in order to maintain con¬ 
centration gradients, but also in order to built and adapt 
cellular structures. Many of these transport issues are 
carried out by molecular motors, i.e. proteins which per¬ 
form a directed stochastic motion along the polar fila¬ 
ments of the cytoskeleton. 

The cytoskeleton of eukaryotic cells is composed of 
three different kinds of filaments: actin, intermediate fil¬ 
aments and microtubules. These filaments give a cell its 
characteristic shape, but also play a key role in intracel¬ 
lular transport. Especially microtubules (MTs), which 
have well-defined plus- and minus-ends, are responsible 
to sustain transport from the nucleus to the periphery 
and vice versa. 

Two families of proteins called molecular motors walk 
along the MTs - kinesin and dynein. The main difference 
between the two is that kinesin walks from the minus- 
end, which is in general located close to the nucleus, to 
the plus-end which grows in the vicinity of the cell mem¬ 
brane, while dynein walks in the opposite direction [5]. 
After this processive motion motors detach stochastically 
from the filament and diffuse around before attaching 
again. 

Molecular motors are able to transport very different 
types of cargo, either individually or by teams of molecu¬ 
lar motors, which are attached to the same cargo. Trans¬ 
port by teams of molecular motors is particularly relevant 
for large objects like vesicles or cell organelles [3]. Trans¬ 
porting cargo by many motors in a crowded environment 
is obviously beneficial since it increases the processivity 
along the MT dramatically and also enhances the ability 
of the motor cargo complex to withstand larger friction 
forces. 

For several types of cargo, however, a bidirectional mo¬ 
tion along a filament was observed in vivo [1] and in vitro 


[ 3 . These observations suggest that the cargo is trans¬ 
ported both by kinesin and dynein motors which are at¬ 
tached to the cargo at the same time. The trajectories of 
these cargos consist of sections of persistent motion and 
sudden returns. One key question is if the observed cargo 
motion is a result of unknown coordination mechanisms 
or is driven by fluctuations. 

Various theoretical models have been suggested, which 
aim at describing the origin of the complex dynamics 
of bidirectionally transported cargo mi- A mean-field 
model (MF-model) which describes bidirectional cargo 
motion [T] driven by two teams of molecular motors, was 
introduced a few years ago by Muller et ai. They as¬ 
sume that two equally strong teams of molecular motors 
with opposed walking directions are bound to a cargo 
at the same time. The model focuses on the force bal¬ 
ance between the two motor-teams without taking the 
motor positions explicitly into account. The forces act¬ 
ing on the motors determine their attachment and de¬ 
tachment rates. The cargo’s velocity depends on the 
forces acting on the motors and is uniquely determined 
for a given number of attached motors of and 
motors. Given the mean-field assumption that forces are 
equally shared among motors moving in the same direc¬ 
tion, the model by Muller et al. predicts the existence of 
high motility states in cargo transport that may originate 
from a purely stochastic tug of war between oppositely 
directed motors, rather than being induced by a regula¬ 
tory mechanism. 

In this work we will test this scenario by using a 
more explicit modeling approach which is inspired by 
the model recently introduced by Kunwar et al. [5]. In 
contrast to the MF-model we explicitly consider the mo¬ 
tor’s positions on the filament, and the couplings between 
motors and cargo which are modeled as linear springs 
in our model (explicit position-based, EPB-model). As 
we want here to test the consequences of the tug-of-war 
mechanism rather than modeling an explicit experimen- 
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tal setup, we consider that and ’’—’’-directed mo¬ 

tors have the same response to applied forces. Our re¬ 
sults show the absence of long-living directed transport 
states. We further extend the model and introduce a 
motor-motor activation, which was inspired by experi¬ 
mental observations [7]. 


Vf 

1000 nm/s 

Fd 

3 

pN 

a 

0.1 pN/nm 

Vb 

6 nm/s 


1 

s-i 

Lo 

110 nm 

Fs 

6 pN 

ka 

5 

s-i 

d 

8 nm 


TABLE I. Simulation parameter taken from [T] and [B]. 


MODEL DESCRIPTION 

In the EPB-model we assume that two teams of motors 
are tightly bound to a cargo. Each team consists of N 
and ’’—’’-motors, respectively. 

To determine the load force applied on the cargo at po¬ 
sition xc(t) at time t we take every single motor position 
Xi into account. We model the motor linker as a linear 
spring with spring constant a and untensioned length 
Lq, such that motors experience no force when located 
at a distance smaller than Lg from xc(t). So the force 
Fi{xc{t), on the cargo caused by the i-th motor is 
given by 


xc{t)) = 


a{Xi - xc{t) + Eg), 

Xi - xc{t) < -Eg 

0, 

\Xi - xc{t)\ < Eg 

a{xi - xcit) - Eg), 

Xi - xc{t) > Eg. 


We assume that no force is exerted if the motor is not 
attached to the filament. When attached the motors can 
perform steps along the filament with a rate depending 
on the force Fi 


where gives the force-free detachment rate and the 
detachment force Fjj gives the force scale. Once a motor 
is detached it attaches again with a rate ka within the 
tensionfree area xc{t) ± Lo¬ 
in the MF-model [1] the cargo moves with a constant 
velocity which is determined only by the number of at¬ 
tached motors of each team between two motor events 
(at-/detachment). To calculate this unique cargo veloc¬ 
ity Muller et al. introduce an artificial force which assures 
a force balance between the two teams of motors with 
equal sharing of the force within one team. All motors 
are implicitly assumed to move with the same constant 
velocity as the cargo. 

In contrast, here we take the individual motor steps 
into account. The force applied on the cargo is deter¬ 
mined separately for each motor (Eq. 0 ) and is modi¬ 
fied when the motor hops. We have thus to define how 
the cargo reacts on the force changes induced by motor 
hopping. We move the cargo of mass m with radius R 
along its equation of motion in a viscous medium with 
viscosity ry [8] 

d^xc{t) ^ jydxc{t) , r 

w—=-67 r?7i?^^^^-h 2^ Fi{xc[t),{x,}). 

(4) 





E, < 0 
0 < Fi < Fs 
Fi > Fs 


( 2 ) 


for -motors and symmetrically for ’’—’’-motors. If 
no force is applied on the motors or if the force is in the 
direction of its motion the motor moves with its force-free 
velocity vp divided by the step length d. ”-|-/—’’-motors 
continue stepping in their preferred direction if the force 
is smaller/bigger than the stall force {+/—)Fs but the 
velocity decreases linearly when the force applied to it 
increases. If the force (in absolute value) exceeds the 
stall force the motors step backwards with vb 'C vp. 

Please note that we don’t consider exclusion of the mo¬ 
tors on the filament, since the number of possible binding 
sites in close proximity of the cargo is huge compared to 
the number of attached motors. 

Next to the stepping rates also the motors’ detachment 
rates kd(Fi) are force-dependent and explicitly given by 

kd{Fi) = k°aexp(^-^'^ , (3) 


Here rip and n_ describe the and ’’—’’-motors, re¬ 
spectively, which are attached to the filament, such that 
0 < < Np and 0 < n_ < N_. 

In the second part of this paper, we shall consider a pos¬ 
sible synchronization of molecular motors via mutual mo¬ 
tor activation in a way we describe later. 


RESULTS 

In this work we analyze whether a model which takes 
every single motor position into account produces the 
same high motility state as it was seen for strong motors 
in the MF-model [1]. Therefore, we consider the same 
kind of strong motors, as defined in TABLE |Tj and we 
measure the velocity distribution and the probability of a 
given number of attached motors of each kind p(n+, n_). 

In the MF-model [1] the cargo velocity is simply a func¬ 
tion of attached motors and is thus constant between two 
attachment/detachment events. In our case the cargo 
moves according to its equation of motion and, there¬ 
fore, not with constant speed. This is why we discretize 
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the trajectory and define v{t) := ^c(t) ^ Here 

we use At = 0.16 s. 


No activation 

In this first subsection we consider no mutual motor ac¬ 
tivation. With the chosen set of parameters given in TA¬ 
BLE]^ the MF-model produces fast cargo motion with a 
bimodal velocity distribution around ±u/, corresponding 
to a bimodal histogram of attached motors (FIG. da)). 
For the EPB-model, however, these high motility states 
are not observed (FIG. db)). This result gives strong 
evidence that fluctuations of the motor positions play a 
crucial role in the cargo dynamics, and that the MF as¬ 
sumption of an equal sharing of forces between motors of 
the same kind qualitatively changes the results. 

For the EPB-model, the distribution of bound motors 
p(n+, n_) in FIG. db) has a peaked structure on the di¬ 
agonal, i.e. most of the time the same number of motors 
of both teams are bound to the filament. Furthermore, 
in the frame of the EPB-model, we can ask how many of 
those bound motors are actually engaged in the tug-of- 
war, i.e. apply a non-zero force to the cargo. This quan¬ 
tity p(n+,n_) shown in FIG.dc) demonstrates that not 
all motors which are bound to the filament exert a force 
to the cargo, a fact which, again, contradicts the mean- 
field assumption. 

When a motor detaches due to a high load, this results 
in the MF-model into a new cargo velocity and a shar¬ 
ing of the high load between the remaining motors of the 
team, which most probably will also detach. In the EPB- 
model, motors have different positions. As illustrated in 
FIG. the motor that detaches is most probably the 
most distant from the cargo, with a reservoir of bound 
motors of the same team behind. After detachment, the 
cargo can move backwards, and as a result, first, the load 
between both teams is slightly relaxed. Second, due to 
the backward motion of the cargo, some bound motors 
which were not exerting any force will become involved 
in the pulling of the cargo and replace the detached mo¬ 
tor. This effect makes obvious why we cannot detect 
the detachment cascades as they are observed within the 
mean-field description, and why high motility states can 
be sustained only in the MF-model. 


Mutual motor activation 



FIG. 1. One specific realization of the trajectories of the 
cargo (blue) and of all (red) and (green) motors. 

The arrows indicate a detachment event and the crosses an 
attachment of a (red) or a (green) motor. The 

color code gives the ratio kd{Fi)/s±{Fi). The initial state is 
chosen such that it would correspond to a high motility state 
in the MF-model, with no ”-|-”-motor attached, and four 
motors pulling the cargo. We see here that this state is not 
stable. 



FIG. 2. The probability p(n+,n_) of n+ ”-|-”-motors and 
n_ ’’—’’-motors bound to the filament (a) in the MF-model 
and (b) for the EPB-model in the case of no activation. The 
probability p(n+, n_) of -|-”-motors and n_ ’’—’’-motors which 
are actually engaged in the tug-of-war (i.e. which apply a non¬ 
zero force on the cargo) in the EPB-model is shown without 
activation in (c) and with a mutual activation with a = 5 and 
Ra ~ 32 nm in (d). 


In the previous discussion we pointed out that, in 
the relevant parameter range, the MF-model can exhibit 
some high motility states (i.e. bimodal distributions of 
the velocity) that are not present in the EPB-model[9]. 
This is because the MF-assumption implicitly introduces 
some synchronization between the motors. We shall now 
explicitely model some mutual motor activation in the 


EPB-model, and explore whether this allows to recover 
a double-peaked histogram of pulling motors as obtained 
in the MF-model. 

The mutual motor activation is introduced as follows. 
If the i-th motor of one team makes a step, the stepping 
rate of the motors of the same team within the interval 
[xi—RA, Xi-\-RA\ is multiplied by a factor a. Here we have 
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(A) (B) 



FIG. 3. The probabilities p(n+, n_) (top) and p(t;) (bottom) 
are shown for a = 1, so without activation on the left hand 
side and for a = 5 and Ra = 32 nm on the right hand side. 

chosen a = 5 and Ra = 32 nm, i.e. 4 tubulin subunits. 

We first analyze the influence of activation for a small 
number TV = 4 of motors of both teams. FIG. id) shows 
that the activation does not change the distribution of 
attached motors significantly. The probability to have 
no motor pulling p(0, 0) goes to zero and those for the 
’’two against one”-scenarios p(2,1) and p(l,2), respec¬ 
tively, are slightly increased. However, we cannot detect 
a double-peaked distribution. 

Now the motor number is increased from N = A io 
TV = 40 motors per team. In column (A) of FIG. 
the probability of engaged motors p(n+,n_) (top) and 
the velocity distribution piv) (bottom) without activa¬ 
tion (a = 1) are shown. The velocity distribution is 
peaked around zero and the maximal probability of en¬ 
gaged motors is still on the diagonal. If we introduce the 
mutual motor activation we observe the bimodal velocity 
distribution and probability of engaged motors as shown 
in column (B) of FIG. Apparently, the mean-field 
assumption is justified for the case of a high number of 
motors and if, in addition, a synchronization mechanism 
between the motors exists. 


DISCUSSION 

In this work we discuss the modeling of motor-induced 
bidirectional cargo transport inside cells. Apart from the 
MF-model of Muller et aT.[T] we also introduce the EPB- 
model which considers explicitly the positions of the mo¬ 
tors bound to the filament in the spirit of [6]. This allows 
us to calculate the forces which act on each individual 
motor explicitly, rather than implicitly via a mean-field 
assumption as in the MF-model. In contrast to the MF- 


model, the EPB-model does not exhibit any high motility 
state corresponding to a bimodal distribution in velocity 
and number of attached motors. 

This result is rather robust and we checked that it does 
not depend on the particular realization of the motor- 
cargo coupling or the details of the motors response to 
external forces. While high motility states are very sta¬ 
ble in the ME-model, a realization with only one team 
pulling the cargo is only a rather rare event with short 
lifetime in the EPB-model. In order to illustrate this 
further we introduced the mutual activation of molecular 
motors as a possible mechanism for synchronization of 
the motors’ motion. When applied to a large number of 
motors, this mechanism led indeed to a stabilization of 
high motility states, comparable to the results of the ME 
model. Therefore, we have shown that fluctuations of 
motor positions are relevant for motor-cargo system, and 
that the mean-field assumption corresponds implicitly to 
an synchronization of the motors which are attached to 
the cargo. A spontaneous synchronization only via the 
tug-of-war between the attached motors can be ruled out. 

In view of this result one has to discuss the experimen¬ 
tal relevance of high motility states. Erom our point of 
view there is no clear evidence that these states have 
actually been observed experimentally, in spite of the 
fact that, for example, the bimodal velocity distributions 
found in have been interpreted as an evidence for high 
motility states. Indeed, in this work, like in many others, 
the trajectories of bidirectionally transported cargos are 
not analyzed directly, but cut into runs and pauses. This 
decomposition of the trajectories follows different dehni- 
tions in different experiments. A common problem of this 
decomposition is that it can induce artifacts in the ve¬ 
locity distribution. Eor example, applying the criteria by 
[To] to the trajectories produced by the EPB-model also 
gives a bimodal velocity distribution (see Supplementary 
Material) though this is not found if we use the full data 
set. More generally, it would be helpful if full trajectory 
data would be systematically provided from experiments. 

Still, several characteristics of cargo transport have 
been found experimentally in in vivo experiments m 
112 ] . including anomalous diffusion (sub- or superdiffu- 
sion depending on the observation time scales). We have 
shown indeed in previous work mE] that, using more 
biologically relevant motor characteristics in the EPB- 
model, these sub- or super diffusive particle motion at 
short times can be reproduced. Still, in vivo experiments 
are not so simple to interpret. In vitro experiments with 
purified motors would allow to have a better control of 
the system parameters, and to go even further in the 
comparison with the models and in the testing of various 
scenarios. 

This work was supported by the Deutsche Porschungs- 
gemeinschaft (DFG) within the collaborative research 
center SFB 1027 and the research training group GRK 
1276. 
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SUPPLEMENTARY MATERIAL 

One quantity to characterize intracellular transport is 
the velocity with which a cargo moves. In experiments 
it is rather difficult to define this velocity. It depends in 
particular on the measurement time window, and there 
is no common procedure to determine it. One popular 
method is to cut trajectories into runs and pauses. Runs 
are usually defined by giving thresholds concerning their 
minimal length, minimal duration and/or minimal veloc¬ 
ity [niin]. Due to the different ways to measure velocity, 
it is rather hard to compare data from different teams. 

Here we want to motivate to use the full data set to 
characterize cargo motion. We have shown in the main 
text that the EPB-model produces a velocity distribution 
which is peaked around zero if we use the whole trajec¬ 
tory to determine it. If we now use the method of m we 
also see a bimodal velocity distribution (FIG.|^. It is ob¬ 
vious why we find this structure: If we ignore the pauses 
the probability to find small velocities goes to zero. So 
by definition the distribution is bimodal. 


From our point of view this kind of determination of 
velocity distributions by only taking parts of the data 
into account can lead to misleading results. 



velocity (nm/s) 

FIG. 4. Velocity distribution of the EPB-model if we use the 
method of [ID] to determine runs and pauses. Here we took 
the parameters of TABLE |I| and A = 4. 

































